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ABSTRACT 

We report on a detailed investigation of the 7-ray emission from 18 broad line radio galaxies (BLRGs) based 
on two years of Fermi Large Area Telescope (LAT) data. We confirm the previously reported detections of 
3C 120 and 3C 111 in the GeV photon energy range; a detailed look at the temporal characteristics of the 
observed 7-ray emission reveals in addition possible flux variability in both sources. No statistically significant 
7-ray detection of the other BLRGs was however found in the considered dataset. Though the sample size 
studied is small, what appears to differentiate 3C 111 and 3C 120 from the BLRGs not yet detected in 7- 
rays is the particularly strong nuclear radio flux. This finding, together with the indications of the 7-ray flux 
variability and a number of other arguments presented, indicate that the GeV emission of BLRGs is most likely 
dominated by the beamed radiation of relativistic jets observed at intermediate viewing angles. In this paper we 
also analyzed a comparison sample of high accretion-rate Seyfert 1 galaxies, which can be considered radio- 
quiet counterparts of BLRGs, and found none were detected in 7-rays. A simple phenomenological hybrid 
model applied for the broad-band emission of the discussed radio-loud and radio-quiet type 1 active galaxies 
suggests that the relative contribution of the nuclear jets to the accreting matter is > 1% on average for BLRGs, 
whilst < 0.1% for Seyfert 1 galaxies. 

Subject headings: Tadiation raechanisms: non-thermal — galaxies: active — galaxies: individual (3C 111, 
3C 120) — galaxies: jets — gamma rays: galaxies — X-rays: galaxies 



1. INTRODUCTION 

A long debated problem in our understanding of accret- 
ing supermassive black holes (SMBHs) in the Universe is the 
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unification of different types of active galactic nuclei (AGN) 
in a framework ascribing their observed diversity to a rela- 
tively few differing parameters and factors. For example, it 
has been widely argued that the difference between type 1 
and type 2 AGN — i.e., those possessing and lacking broad 
permitted emission lines in their nuclear spectra, respectively 
— may be explained by simple geometrical effects involving 
anisotropic obscuration of the active center viewed at different 
inclination with respect to the accretion disk axis ( iAntonuccH 
119931) . The accretion rate was claimed to play a role in this 
context as well, since 'standard' geometrically-thin optically- 
thick accret ion disks formed at high accretion rates (> 1% 
Eddington; Shakura & Sunvaevlll973h on the one hand, and 
radiatively-inefficient geometrically-thick accretion flows ex- 
pected to be present at low accretion rates (< 1% Edding- 
ton; iNaravan & Yil [19941 [T99l lAbramowiczet an[1995h on 
the other hand, can account for different emission spectra and 
power outputs of high- and low-power AGN (se e an early dis- 
cussion on this issue bv'Fabian & Rees'l 995' and'L asotal 19961 
and the more recent one in Ghisellini et al. 2009). The other 
relevant parameter may be also the mass of the SMBH itself, 
which seems to determine some physical differences between 
Narrow-Line Seyfert 1 galaxi es and 'regular' Seyfert Is (e.g., 
iPogge II2OOOI: lKomossall2008l) . 

None of the aforementioned factors can, however, account 
for the dichotomy between radio-loud AGN (which have jets) 
and radio-quiet AGN (which don't). One may therefore 
seek the fundamental difference between AGN producing lu- 
minous radio-emitting outflows and those lacking such out- 
flows in yet another parameter of the accretion disk/black 
hole system. One possibility is that, for example, radio-loud 
AGN harbor rapidly spinning SMBHs with rotational energy 
extracted electromagnetically via the iBlandford & Znaiekl 
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( Il977h mechanism, and converted to the kinetic luminosity 
of relativistic jets. At the same time, the negligible angu- 
lar momentum of SMBHs hosted by radio-quiet AGN pre- 
cludes the formation of such powerful well-collimated out- 
flows. This so-called 'spin paradigm' (Blandford 1990) has 
recently been claimed to be supported, after some mod- 
ifications, by several observational findings and theoreti- 
cal investigations (e . g., iKoide et al.ll2002t ISikora et al.]|2007l: 
lTchekhovskovetal]|2010l see also in this context Garofalo 
2009). Regardless of this debate, the presence of a relativis- 
tic jet constitutes a fundamental distinction between various 
types of AGN, simply because an anisotropic and strongly 
Doppler-boosted jet emission can dramatically affect the ob- 
served properties of a source. In this context, geometrical ef- 
fects play again a major role. In particular, the total radia- 
tive output of those radio-loud AGN which are observed at 
small viewing angles to the jet axis ('blazars') may be domi- 
nated by the broad-band jet emission, while those AGN which 
are inclined at larger viewing angles (e.g., radio galaxies) 
may display radiative signatures of both outflowing a nd ac- 
creting matter at com parable levels (see, e.g., Barthej [19891 : 
lUrry & Padovanil[l99 5). We note that in addition to such ge- 
ometrical effects, the age of a radio structure (i.e., the time 
elapsed after the onset of the jet activity in the nucleus) is yet 
another factor cruci al to understa nding unification of radio- 
loud AGN (see, e.g.. lO'DealfT998h . 

Interestingly, new deep radio surveys indicate that the 
radio-loudness parameter — which is defined as a ratio of 
the jet-related radio flux to the disk-related optical flux, and 
is considered as a proxy for the jet production efficienc\Pl 
— shows a continuous distribution rather than a sharp di- 
vision between radio-loud and radio-quiet AGN. This ap- 
plies to AGN hosted by early-type galaxies and accreting 
at high rates (i.e., quasars), which are typically studied in 
this context (e.g.. White et al. 2000), but holds also when 
other elliptical-hosted AGN (radio galaxies) span ning a wide 
range of the accretion rate are taken into account ( Sikora et alJ 
120071) . Moreover, unresolved non-thermal radio emission and 
jet-like structures have been discovered in classes of AGN 
considered previously as 'radio quiet', i.e., Seyfert galaxies, 
although the jets in such systems are non-relativistic and weak 
when compared to the jets found in 'classical' radio galax- 
ies an d quasars (e. g., Ulvestad & Wilson 1989; Kukul a et alJ 
119951; iThean et aLlJEOOlc iMiddelberg et al.i i2004 and refer- 
ences therein). The distribution of radio-loudness parame- 
ter in Seyferts, which are typically hosted by late-type (disk) 
galaxies is, however, similarly a continuous function of the ac- 
cret i on rat e, as demonstrated first by Ho & Peng (2001) and 
IHoI (I2002h . Yet, it was pointed out by Sikora et^ (i2007i) 
that there is a substantial difference in the distribution of the 
jet production efficiency between disk-hosted and elliptical- 
hosted AGN, with Seyfert galaxies being characterized, at 
any accretion rate, by the radio-loudness parameters orders of 
magnitudes smaller than the analogous parameters character- 
izing elliptical-hosted radio galaxies or quasars. Clearly, more 
studies regarding the jet-disk connection in different types of 
AGN are needed to understand the jet launching processes and 
the physics of active SMBHs in general. 

" More accurately, the radio-loudness parameter depends not only on the 
total power of a jet relative to the total power of an accreting matter, but also 
on their respective radiative efficiencies. Only if these radiative efficiencies 
are not very different between various types of AGN (and various AGN of the 
same type), the radio-loudness parameter may be considered as a good proxy 
for the jet production efficiency. 



In this context, broad line radio galaxies (BLRGs) seem 
ideal targets for an in-depth investigation, since this particu- 
lar class of very radio-loud AGN exhibits both the disk-related 
('Seyfert-like') and the jet-related ('blazar-like') radiative sig- 
natures in their broad-band spectra. Unlike blazars, the jets 
in BLRGs are not pointing directly toward the observer, and 
so the relativistic beaming effects and the related jet domi- 
nance are only moderate. Moreover, unlike narrow-line radio 
galaxies (NRLGs) — which are believed to be intrinsically 
similar but simply inclined at systematically larger jet view- 
ing angles — BLRGs are not generally obscured by large 
amounts of dust distributed in torus-like structures around 
the nucleus, and hence radiative properties of the accretion 
disks and of the circumnuclear gas can be easily accessed 
in their case. BLRGs show in particular optical/UV contin- 
uum and emission-line characteristics very similar to those 
of luminous Seyfert galaxies, which indicates high accretion 
rates and the presence of standard Shakura-Sunyaev disks 
in both classes of objects. Some fundamental differences 
in the X-ray spectra between BLRGs and high-accretion- 
rate Seyferts have been however noted. Specifically, even 
though the observed X-ray/soft 7-ray emission of BLRGs 
seem still dominated by the moderately absorbed emission 
by the accreting plasma (i.e., disks and disk coronae), rather 
than by the jets, the 1 — 100 keV continua of BLRGs are 
flatter, and their reflection components (as well as the flu- 
orescent Fe Ka lines) weaker than in the case of luminouj^ 
Seyfert galaxies (e.g ., Maraschi et al. 1991; Wozniaket 
1998; Sambruna et al. 1999, 2002; Eracleous et alj 12001 
Zdziarski & Grandi 2001 ; Ballantyne 2007 ). Because of 
such differences, several authors in the past speculated about 
the non-negligible jet contribution to the X-ray emission of 
BLRGs, diluting the accretion-related radiative output in the 
X-ray domain (Woz niaket aP 119981: lEracleous et aP 120001: 
Grandi et al. 2002). This idea was subsequently examined 
in various different approaches using most recent broad- 
band X-ray data obtained with BeppoSAX dGrandi et al.ll2006l: 
l&andi & Palumbo 2007), Suzaku and Swift dKataoka et al] 
7007; Sambruna et al. 200|l. 

Grandi & Palumbo (200^ attempted to disentangle the jet 
and the disk contributions to the X-ray spectra of three of the 
brightest BLRGs (t he approach first applied to the case of the 
quasar 3C 273; see iGrandi & Palum bo 2004). For simplicity, 
they assumed that the accretion disks in BLRGs and Seyfert 
1 galaxies produce similar emission continua and reprocess- 
ing features, and subsequently allowed for a presence of the 
Doppler-enhanced jet radiation at an arbitrary level. The fits 
obtained for 3C 120, 3C 390.3 and 3C 382 showed that the 
data are indeed consistent with a combination of a thermal 
component (in a first approximation associated with an ac- 
cretion disk) and a non-thermal component associated with 
the beamed radiation (due to a jet). Grandi & Palumbo con- 
cluded however that jets make only minimal contribution to 
the X-ray continuum e mission of BLR Gs, in agreement with 
the previous findings b y lWozniak et akl (IT998). More recently, 
ISambruna et al.l (12009!) also proposed that BLRGs may be just 
clustered at one end of the distribution of the X-ray spec- 
tral parameters (e.g., photon indices and reflection albedos) 
characterizing Seyfert galaxies with significant overlap. Inter- 
estingly, both Grandi & Palumbo (2007) and Sambrun a et alJ 
02009.) suggested that the emission of the underlying jet may 
instead dominate the radiative output of BLRGs at high en- 
ergy 7 rays, and in particular in the GeV regime. 

With the successful launch of the Fermi Gamma-ray Space 
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Telescope, we now have an unprecedented opportunity to 
study in detail the 7-ray emission from different types of 
extragala ctic so urces — not only b l azars, but also radi o 
galaxies (lAbdo" et al. 20091IH l2010clldl: iKataoka et alllMoh . 
and other classes of AGN as well (such as Narrow-Line 
Seyfert 1 galaxies, for example: lAbdo et al. 2009d) . Dur- 
ing the first 15-month of the Fermi mission, 11 non-blazar 
AGN have been de tected in the GeV photon energy range 
JAbdoetaLl l2010alldl) bv t he Fermi Large A rea Telescope 
rLAT: lAtwood et alJ 120091: lAbdo et alJ l20T0 b). This 'mis- 
aligned AGN sample' includes seven Faranoff-Riley type I 
(low-power) radio galaxies, and four Faranoff-Riley type II 
(high-power, hereafter 'PR IF) radio sources consisting of 
two radio galaxies and two steep spectrum radio quasars. 
The two luminous radio galaxies detected in 7 rays are, 
in fact. X-ray bright objects classified spectroscopically as 
BLRGs: 3C 120 (FR I radio morphology), detected for the 
first time with Ferm/-LAT , and 3C 111 (FR II radio morphol- 
ogy), whose 7-ray de tection was initially reported by EGRET 
jHartman et al.ll2008h and confirmed by Ferm/-LAT . None of 
the X-ray bright Seyfert 1 galaxies appear to be detected in 7 
rays thus far. 

In this paper we report on a detailed investigation of the 7- 
ray emission from 18 hard X-ray brightest BLRGs, as well as 
a comparison sample of high accretion-rate Seyfert 1 galax- 
ies selected as their supposed radio-quiet counterparts (in the 
framework of the AGN unification scheme). Our primary 
goals are to examine the 7-ray properties of BLRGs as po- 
tential '7-ray loud' active galactic nuclei, to study the high- 
energy jet emission in the selected sources, in particular, in- 
vestigating the relative contributions of the nuclear jet and ac- 
cretion disk emission in the broad-band spectra of BLRGs. 
The two years of ferm/-LAT exposure provides us with a 
rather deep all-sky survey reaching the flux limit of typi- 
cally a few xlO^^^ergcm^^ s^^ (95% confidence level) be- 
tween the observed photon energies lOOMeV and 10 GeV. 
In §2, we describe the FermZ-LAT observations and data 
reduction procedure. The results of the analysis are given 
in §3. The discussion and final conclusions are presented 
in §4. Throughout this paper, a ACDM cosmology with 
Ho = 7 1kms~^Mpc-\ O a = 0.73, and n„, = 0.27 is 
adopted JKomatsu et alll2009l) . 

2. DATA AND DATA ANALYSIS 

2.1. The Sample 

Our sample includes all the BLRGs observed by modern 
X-ray astronomy satellites (EXOSAT, Ginga, ASCA, RXTE, 
BeppoSAX, Chandra, XMM-Newton, Suzaku, INTEGRAL and 
Swift), for which data are available at energies above 2 keV. 
Table 1 presents the list of 1 8 BLRGs compiled, which also 
includes for comparison a sample of ni ne bright Seyfert 1 
gal axies chosen from the compilation by iHo & Pen3 (1200 lb 
and lHol (120021) . for which the nuclear optical and hard X-ray 
fluxes match those of the discussed BLRGs. Here we have 
selected luminous Seyferts 1 with measured radio fluxes (or 
meaningful upper limits) which, due to their high accretion 
rates (> 1% Eddington) and unobscured nuclei, may be con- 
sidered as radio-quiet ana logues of BLRGs (see the discus- 
sion in iSikora et"an I2007IP1 . The main difference between 
the two analyzed classes of sources that should be empha- 

^^ In our sample of Seyfert galaxies we have therefore excluded type 1.5 — 
2 objects, in order to avoid additional complications related to the absorption 
of the X-ray emission by the circumnuclear cold gas and dust. We have also 



sized once more is that Seyferts are radio-quiet and hosted 
by late-type galaxies, while BLRGs are very radio-loud and 
elliptical-hosted. The other possibly relevant (but related 
to the morphologies of host galaxies) difference regards the 
masses of their SMBHs, A^bh- In particular, nine Seyferts 
included in our sample are characterized by lower values of 
black hole masses (median « IO^-^AIq) when compared to 
the ten BLRGs considered here with A^bh provided in the 
literature (median « IQ^-^Mq; see Table 1). 

The basic information collected from the literature regard- 
ing each analyzed source, as listed in Table 1, are (1) lAU 
coordinates for J2000, (2) source name, (3) redshift z, (4) 
luminosity distance c^l, (5) total 5 GHz flux [j^fiy]5°GHz ^^ 
the cgs units of ergcm^^s^^, (6) 5 GHz flux of the unre- 
solved nucleus [j^^iylscfHz' C^) -S-band optical flux of the nu- 
cleus [i/F^]g"'^, (8) X-ray photon index Fx measured between 
2 and lOkeV, (9) average 2 - lOkeV flux [i^F^]2-iokcV, 
(10) hard X-ray/soft 7-ray flux [i^Fy]i4-i95kcV detected by 
Swift-BAT at 14 - 195 keV, (11) black hole mass Mbh, and 
(12) references. In the case of Seyf erts the provi ded nu- 
clear S-band fluxes (same as in Sikor a et aP l2007h . which 
are carefully corrected fo r the non -negligible starlight con- 
tamination, are taken from lHo & Pe ng (2001) and Ho (2002). 
In the case of BLRGs, much less severe starlight contami- 
nation was taken into accoun t by means of the appropriate 
correction factors give n by Er acleous & Halpern I (1994) and 
lEracleous & Halpern 1 020031) . Also, when a _B-band flux was 
not provided in the literature explicitly, we estimated it from a 
V-band flux applying a multiplication factor (Av/Ab)^"""^' 
with Ay = 5500 A, Ab = 4400 A, and the assumed optical 
power-law slope aopt = 0.5. Finally, all the S-band nu- 
clear fluxes listed in column 7 have been corrected for the 
Galactic extinction available in NEEE^. We note that BLRG 
3C 227 has been detected only very r ecently by Sw ift-BAI at 
a high significance level (S/N of 4.87: lTueller et ay2010 ). but 
is poorly known at other wavelength. In X-rays, two Chan- 
dra observations of this object have been conducted so far, as 
parts of a survey of multiple hot spots in the large-scale struc- 
tures of nearby radio galaxies. We estimated the 2 — lOkeV 
nuclear flux of 3C 227 from the count rate of Chandra ACIS 
CCD chip given in lHardcastle et al] (l2007h using the software 
PIMMS0 and assuming the X-ray photon index Fx = 15. 


We are aware that the compiled sample of BLRGs is by 
no means statistically complete or unbiased, since the tar- 
gets were originally selected from independent programs by 
different observers using different instruments. The final 
list of objects reflects in fact a bias toward sources which 
are the brightest in the X-ray domain. Moreover, we are 
aware that most of the considered objects s how flux variabil- 
ity in various energy bands. For example, iTurner & Poundsl 
(Il989h recorded the minimum X-ray flux of 3C 111 at the 
level of 1.8 x lO^^^ergcm^'^s^^, whereas the maximum 
flux for this source recorded a decade later reached 5.6 x 
10~^ ^ ergcm^-^ s^^ in the 2 — lOkeV band (Eracleo us et alj 
l2000h . Table 1 lists the average fluxes measured at different 



exclude d sources such as N GC 4639, which, even though classified as type 1 
Seyfert fSikora et al."200T), is an example of a low-luminosity AGN accreting 
at a low rate (Ho et al. 1999). 

http : //nedwww. ipac . caltech . edu 
http ://heasarc.nasa. gov/ Tool s/w3pimm. html 
^' The uncertainty in the assumed X-ray photon index ±0.5 would make 
only ~ 20% difference in the estimated 2 — 10 keV flux. 
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Fig. 1 .— Temporal variation of 7-ray flux (> 100 MeV) of 3C 111 (top) 
and 3C 120 (bottom) over the period 2008 August-2010 August. The time 
(in days) is measured from the start of the Fermi'-LAT observation (2008 Au- 
gust 4, 15:43:38 UT). The fluxes are plotted only when the coiresponding TS 
values exceed 10 for a given time bin, otherwise upper limits are provided. 



wavelengths as provided in literature, while at the same time 
the results discussed later in this paper are not affected by 
the relatively modest (typically up to a factor of 3) flux vari- 
ations characterizing X-ray and radio continua of BLRGJ^. 
The same caveats apply to the sample of Seyfert galaxies ana- 
lyzed as it was compiled using even more heterogeneous cri- 
teria. Hence it is difficult to identify and to discuss all the 
possible biases introduced by the applied selection criteria for 
both samples. Therefore these should be considered simply 
as lists of prominent (bright) examples of the two discussed 
classes of AGN. 

2.2. Fermi-LAT Observations and Data Analysis 
The LAT instrument onboard Fermi is described in detail in 



lAtwood eT al. (2009) and references therein. It is character- 
ized by a larger effective area (~8,000cm^ on axis at 1 GeV 
for the event class considered here), a wider energy coverage 
(from ^20 MeV to >300GeV), and an improved angular res- 
olution when compared to the previous 7-ray missions. The 

-'' In order to illustrate possible effects caused by the temporal variability 
of the analyzed sources, in the spectral energy distributions shown in Fig- 
ures 3, 8 and 9 below we have provided both the minimum and maximum 
X-ray fluxes whenever available in the literature. 



68% containment angles of the reconstructed incoming pho- 
ton direction are approximated as Sgs — 0°.8 (e^/GeV)^"'® 
below 10 GeV. During the first two years of FermZ-LAT oper- 
ation, most of the telescope's time was dedicated to observing 
in a 'survey mode,' in which the instrument points away from 
the Earth and nominally rocks the spacecraft axis north and 
south from the orbital plane to enable monitoring of the entire 
sky on a time scale shorter than a day. In particular, the whole 
sky is surveyed every ~ 3 hours (2 orbits). 

The dataset used here comprises all the scientific data ob- 
tained between August 4, 2008 and August 4, 2010. This time 
interval runs from Mission Elapsed Time (MET) 239557414 
to 302630530, which is consistent with the observation period 
for the Second Fermi-hPC[ Catalog (2FGL) selection (The 
Fermi-LAT collaboration 201 1, in prep). We have applied the 
zenith angle cut of 105° to eliminate photons from the Earth's 
limb. We use events from the "Diffuse" class (Atwoo det al] 
2009), i.e., the events that have the high probability of be- 
ing photons. In the analysis presented here, we set the lower 
and higher energy bounds at 200 MeV and 100 GeV, respec- 
tively. The choice of a lower energy bound at 200 MeV is 
conservative but reduces significantly systematic errors. Sci- 
ence Tool0 version v9r15p2 and Instrumental Response 
Functions (IRFs) P6_V3_DIFFUSE were used throughout 
the analysis. 

In order to study the average spectrum of each selected tar- 
get, we use the standard unbinned maximum-likelihood spec- 
tral estimator (Mattox et al. 1996) provided with the LAT sci- 
ence tool GTLIKE. It allows US to fit the data to a source 
model, along with the models for the uniform extragalactic 
and structured Galactic background^ The Galactic diffuse 
emission model and the isotropic spectral model used here 
were developed by the Ferm/-LAT team as refinements of the 
publicly released models; this choice does not affect signifi- 
cantly the results for the candidate sources considered, which 
are all located outside the Galactic plane except for 3C 1 1 lE3- 
Much more crucial is a careful selection of source regions, 
especially for relatively faint objects. The model for which 
we calculate the likelihood is a combination of point-like and 
diffuse sources for a region of interest (ROI) with the radius 
of r = 8° centered on the target under consideratioro For 
all the BLRGs and Seyferts listed in Table 1, we first as- 
sume a point-like source with a power-law spectrum at the 
position of each target, and fix the 7-ray photon index as 
r^ = 2.5. Additional sources f rom the 1 st Ferm/-LAT Cat- 
alog (IFGL: lAbdo et all 1201 Obi) and from the internal LAT 
collaboration catalog produced using 18 months of data are 
included in the model of each ROI. Next, using the GTLIKE 
tool, we find the best-fit parameters for each source and evalu- 
ate the significance of the detection given by the test statistic, 
TS = 2A log(likclihood) between the models with and with- 
out a source. For sources which are detected above a certain 
significance threshold (TS > 25, corresponding to > ba de- 
tections), we retry the GTLIKE fit assuming power-law spectra 
(F = K E^^'') with both parameters (normalization, K and 

http: //fermi . gsf c . nasa . gov/ssc/data/analysis/documentation 
http ; //fermi . gsf c . nasa . gov/ssc/data/access/lat/BackgroundM 
2' With Galactic coordinates (I = 161°. 68, b = -8°. 82) 3C 111 is 
located behind a large molecular cloud complex in Taurus (cf. lDame et alj 
1200 IL Fig. 2 therein), thus requires more care in the analysis due to the ex- 
pected contamination from Galactic diffuse emission. 

^* See, in this context, Abdo et al. (2009a) for a detailed investigation of 
changing the radius of ROI from 8° to 20° , and the arguments given for using 
r = 8° to minimize the contamination from the Galactic diffuse emission. 
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Fig. 2. — Fermi-LAT TS map (>200 MeV) centered on Pictor A, show- 
ing the presence of multiple 7-ray peaks in the field. Contamination from 
nearby sources ai'e modeled as part of the background that also includes 
Galactic/extragalactic diffuse 7-ray emission. The green contours corre- 
spond to TS values of 5, 10, 15. and 20. while the red contours indicate 1.4 
GHz radio emission from Pictor A iPerley et al. 1997). The position of a 
blazar, BZQJ0515-4556, is also shown as red box. The peak near the center 
of the map (TS = 20) is almost exactly coincident with the position of Pictor 
A. The TS peaks associated with Pictor A and with BZQ J0515— 4556 ai'e 
only marginally resolved. The TS value quoted here and in Table 2 for Pic- 
tor A was evaluated using a source model that included point sources at the 
position of BZQ J0515-4556 and at R.A. = 81.35 deg, Dec. = -45.78 deg. 

photon index, Tj) set free, and then calculate the errors on 
the fluxes and photon indices. For the remaining BLRGs and 
Seyferts, which are detected in the analyzed dataset below the 
threshold (TS < 25), we simply provide upper limits on the 
fluxes for the fixed F^ = 2.5. 

3. RESULTS 

Table 2 summarizes the results of the ferwi-LAT data anal- 
ysis of the 18 BLRGs and 9 Seyfert galaxies. For each source 
considered. Table 2 provides (1) name, (2) statistical signif- 
icance TS of the Ferm/-LAT detection, (3) 7-ray photon in- 
dex F-y evaluated for the photon energy range 0.1 — lOGeV, 
(4) the integrated photon flux above 100 MeV, F>o.igcV. (5) 
7-ray flux [i/Fi,]o.i-io GeV> (6) the corresponding 7-ray lumi- 
nosity Lj = A7rd^[i'F,y]o,i-ioGcV7 (7) total accretion-related 
luminosity Lace derived from the spectral fitting as described 
below, and (8) 'mixing' parameter 77 discussed in the next sec- 
tion. Only two BLRGs (3C 111 and 3C 120) are detected at 
sufficiently high significance levels, i.e., TS > 25, in the ac- 
cumulated two-year Ferm/-LAT dataset. For these, the 7-ray 
fluxes and luminosities are evaluated straightforwardljEj- For 
the targets detected at lower significance levels, TS < 25, 
the corresponding 95% confidence level flux upper limits are 
calculated using the dedicated software UpperLimit.PY. 

For 3C 120, th e results presented h ere are consistent with 
those reported in lAbdo et aP (l2010dl) . but with reduced un- 
certainties in the flux and photon index due to the improved 
photon statistics based on the two-year accumulation of the 
Fen«/-LAT data. We note however that the TS value in- 
creased only slightly between the 15 -month and 24-month 
Fermi-L AT datasets (cf., TS = 34 found here vs. TS = 32 re- 
ported in lAbdo et al.ll2010d '). In contrast, the flux and photon 
index uncertainties increased in the case of 3C 111, and the 
corresponding TS value decreased between the 15-month and 

^' Since the likelihood analysis was limited to the photon energy range 
0.2 — 100 GeV, all the flux values and the coiTesponding luminosities are 
extrapolated down to 100 MeV with a given photon index. This choice is 
dictated solely by the convention typically followed by the 7-ray community. 
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Fig. 3. — Broad-band SEDs of the two BLRGs detected at high significance 
by Fermi-hAI (3C 111 and 3C 120). Fermi-hAI data are indicated by red cir- 
cles. Black squares represent the historical data from NED. Magenta squares 
denote the 5 GHz radio fluxes of the um'esolved nuclei (if available). The 
green curves correspond to the template of the accretion-related Seyfert-type 
emission (from Koratkar & Blaes 1999), matched to the infrared-to-X-ray 
continuum of each source. The blue curves correspond to the broad-band 
spectrum of the quasar 3C 273 (from Soldi et al. 2008), used here as a tem- 
plate of the jet-related emission and scaled to match the radio fluxes for each 
source. The mixing parameter r] for the phenomenological hybrid model dis- 
cussed in § 4 is given in each panel. 



24-month Ferm/-LA T datasets (cf. TS = 31 found here vs. 
TS = 34 reported in lAbdoet al.ll2010dl) . The reason for this 
behavior is twofold. First, the likelihood analysis was lim- 
ited here to the photon energy range 0.2 — 100 G eV, whereas 
the ene rgy range 0.1 — 100 GeV was adopted in lAbdo et alj 
(!2010d). The difference in energy selection is relevant since 
3C 111 is located at a relatively low Galactic latitude (see 
footnote |27| i, and as such is heavily affected by the contam- 
ination from the Galactic diffuse emission especially below 
200 MeV H If we lower the photon energy threshold of the 
likelihood analysis to 100 MeV, the significance of the 3C 1 1 1 

'" Accordingly, in the IFGL catalog I Abdo et al. 2010hl table 4 therein), 
the source fit of IFGL J0419. 0-1-3811 was flagged as being sensitive to 
changes in the diffuse Galactic emission model (flux and spectral index could 
change by more than 3a-). However, upon close inspection of the dust col- 
umn d ensity E{B~V) and ly(CO) maps I Schle gel et al.1199 81: Dame et al.l 
1200 Ih . we estimated that the possible enhancement of 7-ray emission at the 
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Fig. 4. — Distribution of a mixing parameter rj for BLRGs {top panel; blue 
histogram) and Seyferts (bottom panel; red histogram). The values charac- 
terizing the two detected sources (3C 111 and 3C 120) are indicated by the 
arrows in the top panel. In the case of Fairall 9, which is not detected in radio, 
an upper limit for the rj parameter is given in the bottom panel. 



detection in the 24-month d ata increases (TS = 59), as ex- 
pected. Second, as noted in ' Abdo et al] (l2010dl) . the consid- 
ered galaxy seemed variable in the GeV range, being in par- 
ticular brighter at the very beginning of Ferm/-LAT observa- 
tions. Here we investigate this issue in more detail, showing 
in Figure 1 (top panel) the temporal variations of the 7-ray 
flux of 3C 111 above 100 MeV using the two-year accumu- 
lation of the Ferm/-LAT data binned in the 3-month-long pe- 
riods. The fluxes are plotted only when the detection signif- 
icance reached TS > 10 in a given time bin; in the case of 
3C 1 1 1 such a criterion was fulfilled only in three time in- 
tervals out of 8 total. Even lower 7-ray duty cycle emerges 
from the 3C 120 lightcurve, as also shown in Figure 1 (bot- 
tom panel). All in all, we conclude that both BLRGs detected 
by ferm/-LAT are variable in the GeV photon energy range, 
even though the significance of the variability can hardly be 
evaluated due to the insufficient photon statistics. A dedicated 
analysis of all EGRET data also indicate d plausible variabilit y 
of the MeV/GeV emission in 3C 111 (Hartman et al.ll200a) . 
Note also that over the past two years both BLRGs have de- 
clined in flux by ^30% at centimeter wavelengths as observed 
by the Univ. of Michigan Radio Astronomy Observator}0- 
In 3C 111 specifically, a bright sub-mm (230 GHz) flare was 
observed toward the end of 2008 ( Chatteriee et alj|201 Ih . co- 
inciding with the period of the initial ~0.5 years of LAT ob- 
servations, with a subsequent decline in the sub-mm flux over 
the next ^1.5 years. 

At this point let us also comment on the particular case 
of Pictor A galaxy, for which a relatively high TS value has 
been found (see Table 2), yet below the critical value of 25 
required for claiming the detection. In fact, in the analysis 
procedure described above (i.e., assuming a single point-like 
source at the position of the galaxy, R.A. = 79°. 957 and 
Dec = -45°. 779), one can find formally TS = 45 for Pic- 
tor A. However, a close inspection of the spatial distribution 
of TS for the whole source region (hereafter the 'TS map') 
reveals a rather complex structure elongated substantially in 
the RA direction, which is inconsistent with a single point 

position of 3C 111 due to additional column density from dark gas or a some- 
what larger emissivity may account for at most, ^ 10 % of the 7-ray flux listed 
in Table 2. 

http://www.astro.lsa.umich.edu/obs/radiotel/umrao.php 
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Fig. 5. — Top panel: the optical (B-band) versus X-ray (2 — lOkeV) fluxes 
for all the sources included in the sample, with BLRGs denoted by blue filled 
circles and Seyferts by red crosses. Bottom panel: the radio (5 GHz) versus 
X-ray (2 ~ lOkeV) fluxes. Here BLRGs are plotted as blue filled cycles 
when nuclear radio fluxes are used, and as cyan open circles when total radio 
fluxes are considered. Similarly, Seyfert 1 galaxies are plotted as red crosses 
when nuclear radio fluxes are used, and as magenta crosses when total radio 
fluxes are considered. An upper limit is given for Fairall 9. In both panels, 
large symbols indicate the two objects detected by Fermi-hAI (3C 111 and 
3C 120), while the medium-size symbols represent Pictor A galaxy which, 
even though not formally detected, is characterized by the third-highest TS in 
the analyzed two-year fenni-LAT dataset. 



source. Figure 2 shows TS map thus obtained, where contam- 
ination from nearby sources (listed in the IFGL and in the in- 
ternal LAT collaboration catalog using 18 months of data) are 
modeled as part of the background that also includes Galac- 
tic/extragalactic diffuse 7-ray emission. This suggests a pres- 
ence of multiple 7-ray sources in the field. More exactly, in 
a careful examination of the TS map we found three emission 
peaks (each about a degree apart), one of which coincides al- 
most exactly with the position of Pictor A, being character- 
ized by TS ~ 20. Also, one of the other two emission peaks 
located at R.A. = 79°. 09, Dec = -46°. 08 is positionafly 
coincident with a blazar BZQ J0515-4556 (R.A. = 78°.94, 
Dec — —45°. 95) with a TS value of 19. The last peak seen 
to the east of Pictor A is less prominent (TS value of 1 1) and 
located at R.A. ^ 81°. 35, Dec = -45°. 78. Sources at the 
positions of the two TS peaks just mentioned were included 
in the model when the TS for Pictor A listed in Table 2 was 
evaluated. Hence, one may conclude that there are some indi- 
cations for the GeV emission of Pictor A galaxy close to the 
upper limits provided in Table 2, and that the formal detection 
of this object by Fermi-LAI in the near future is quite likely. 
No statistically significant detection of any Seyfert 1 galaxy 
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analyzed was found in the 24-month Fermi-LACT dataset. In- 
terestingly, IC 4329 A shows a relatively high TS value of 15 
(see Table 2), but this could be due to a contamination from a 
nearby source. This nearby source was found only recently, 
and has been associated tentatively with a blazar in the 2FGL 
catalog (The Fermi-LAT collaboration 2011, in preparation). 
Nevertheless, the provided upper limits for the 7-ray emission 
of the considered sample of Seyferts — typically at the flux 
level of a few xlO^^^ergcm^^ s^^ between 0.1 and lOGeV, 
depending on the degree of the contamination from the Galac- 
tic diffuse emission and on the presence of nearby bright 7-ray 
sources — provide interesting constraints as discussed below. 

4. DISCUSSION AND CONCLUSIONS 

The analysis of the two year Fermi-LAT data for the se- 
lected 18 X-ray-bright BLRGs confirmed the previously re- 
ported detections of the two sources (3C 120 and 3C 111), 
at the significance levels not much different between the 
15-month and 24-month datasets. This may indicate that 
the observed 7-ray emission is variable on months to year 
timescales, with a relatively low duty cycle, as revealed by a 
closer inspection of the corresponding lightcurves (Figure 1). 
In fact, the > 100 MeV 7-ray flux of 3C 1 1 1 observed with 
Fermi-LAT is (3.5±1.2)x 10^^ ph cm^^ s^^, which is about 
~ 20x smaller than the maximum recorded by EGRET in the 
same energy range (64x10^^ ph cm^^ s^^, iHartman et al.l 
I2OO8I) . thus suggesting significant variability in the decade 
timescale between the EGRET and Fermi-LAT observations. 
Moreover, we found some hints for 7-ray emission from Pic- 
tor A (TS — 20), even though we cannot claim a formal 
detection at the moment. These results suggest that other 
BLRGs could be promising candidates for Fermi-LAJi detec- 
tions in the near future, and that BLRGs in general could po- 
tentially constitute a class of 7-ray loud AGN. Such a state- 
ment is however difficult to quantify, because of the afore- 
mentioned variability of the GeV fluxes from 3C 120 and 
3C 1 1 1 (assuming that both galaxies are representative for the 
whole class). On the other hand, the observed flux changes 
from these two detected galaxies imply that, in analogy with 
blazar sources, the GeV continua of BLRGs are produced pre- 
dominantly in the inner parts of their nuclear outflows (jets 
on scales less than tens or hundreds of parsecs). In fact, as 



noted in the previous sections, several authors have previ- 
ously suggested that the GeV radiation from BLRGs, if de- 
tected, should most likely be due to non-thermal and beamed 
blazar-type emission, but only observed at intermediate incli- 
nations (jet viewing angles, say, 10° < 64 < 30°, versus 
Oj 5, 10° expected fo r blazars; see Grandi & PalumboluOO Tc 
ISambruna et al.ll2009l) . Interestingly, unlike in the luminous 
blazars, the X-ray/soft 7-ray spectra of BLRGs (up to the ob- 
served photon energies of ^ 100 keV) have been argued to 
be dominated by the thermal radiation of the accreting matter 
([Wozniak e t al. 1998; Sambruna et al. 1999; Eracleous et al] 
I2OOOI: iZdziar ski & Grandi 200lJ) with only little (if any) jet 
contribution CGrandi & Palumbol lIOOTt iKataoka et al.1 120071: 
ISambruna et al.ll2009h . Hence it is clear that BLRGs are truly 
ideal targets for investigating the AGN jet-disk connection in 
the X-ray/7-ray regime. Here we present a few considerations 
regarding this issue. 

Figure 3 presents spectral energy distributions (SEDs) of 
the two aforementioned 7-ray-detected BLRGs; the SEDs 
of all the remaining objects, for which only the upper lim- 
its in the Fermi-LAil range were derived, are given in Fig- 
ures 8 and 9 further below (see Appendix A). In the figures, 
the Ferm/-LAT data are indicated by red circles (or arrows), 
black squares represent the historical data from NED, while 
the magenta squares show the 5 GHz radio fluxes of the un- 
resolved nuclei. Even though the broad-band spectra vary to 
some degree between 3C 120 and 3C 1 1 1 (e.g., with respect to 
the ratio between the total and the nuclear radio luminosities), 
one can gauge the main similarities and differences between 
the BLRG-type and Seyfert-type SEDs. In particular, while 
the infrared-to-hard X-ray continua of the two BLRGs are 
very similar to the one characterizing Seyfert galaxies consid- 
ered here, the former objects are significantly brighter than the 
representative Seyfert in the radio and 7-ray domains. Such 
a dramatic difference in the radio loudness is related to the 
presence of a relativistic jet, as discussed above, and here we 
suggest that same is true regarding the 7-ray loudness. To in- 
vestigate it further more quantitatively (yet still illustrative), 
we consider a simple phenomenological 'hybrid' model for 
the broad-band emission of a type 1 AGN, consisting of indi- 
vidual thermal and the non-thermal emission components (the 
approach analogous to the one adopted in Grandi & Palumbol 
2004, 2007). The thermal component, related to the radiative 
output of the accreting and circumnuclear matter (accretion 
disk, disk corona, dusty torus), is approximated here by th e 
template Seyfert spectrum given bv lKoratkar & BlaesI (II999I) . 
For the non-thermal broad-band emission component of the 
nuclear (blazar-type) relativistic jet observed at intermediate 
viewing angles, we use the well-constrained SED of the radio- 
loud quasar 3C 273 (from lSoldiet al.ll2008h as a templateQ 
Both templates, plotted in Figure 3 (as well as in Figures 8 
and 9) as green and blue curves, respectively, are rescaled to 
match the data points for the analyzed sources. If the Fermi- 

^^ It should be emphasized here that the observed multiwavelength 
emission of 3C 273 is not purely non-thermal in origin. For example, 
Grandi & Palumbo (2004) estimated that about a quarter to half of the 
2 — 10 keV flux of this source could be due to the accreting matter, and not 
due to the jet. However, the "big blue bump" is clearly visible in the overall 
SED of the quasar around UV/soft X-ray frequencies is almost certainly due 
to the accretion disk. Nevertheless, the broad-band spectrum of 3C 273 is 
the best sampled (and best understood) spectrum of a quasar dominated in 
radio and in 7 rays by the moderately-beamed emission of a powerful jet ob- 
served at intermediate viewing angles (9j ~ 10°; see Courvoisiet 1998, and 
references therein). As such, it can indeed be consider as the best available 
template for the jet-related emission of BLRGs. 
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BLRGs detected by Fermi-LAJ (3C 111 and 3C 120), while the medium- 
size open circles represent Pictor A. In the bottom left panel, when total radio 



LAT did not detect the source and only upper limits on 7-ray 
flux were derived, its non-thermal component is entirely de- 
termined by the 5 GHz nuclear flux (in this context, see Ta- 
ble 3 of Ghisellini et al. (2005), who predicted 7-ray fluxes of 
3CR FR I radio galaxies from 5 GHz core fluxes). This al- 
lows us to approximate the relative contribution of the non- 
thermal and thermal emission components to the observed 
SED of each object, in terms of a 'mixing' parameter (77). Un- 
der the adopted definition, this parameter increases with the 
increasing contribution of the jet-like component, and equals 
unity when both the 3C 273 and the Seyfert template spec- 
tra provide comparable contributions to the observed UV flux 
of a source. Note that the direct radiation of the standard 
(Shakura-Sunyaev) AGN accretion disk is expected to peak 
in the UV regime (photon energies ^ 10 eV). 

Obviously, the adopted model is an oversimplification of a 
realistic situation, as it ignores several complications which 
may be potentially relevant. For example, it is at some level 
questionable to use the broad-band spectrum of the radio- 
loud quasar 3C 273 as a template for the jet-related emission 
of Seyfert galaxies in general. Moreover, possible (or even 
likely) temporal variability of each source considered, as well 
as the ad hoc adopted procedure of matching the assumed 
templates to the collected data points, precludes any robust 
statistical analysis (e.g., x^ fitting) which would allow for the 
values and errors of the 77 parameter to be determined more 
accurately. Nevertheless, it is worth noting that this simple 
model seems to match the SEDs of the two BLRGs well from 
radio to 7-ray frequencies, and that in both cases the emerging 
values of the 77 parameter are similar (~ 0.15 — 0.35), being in 
ad dition consistent with the o nes following from the analysis 
bv iGrandi & Palumbol (l2007h . Importantly, the model when 
applied to the other BLRGs analyzed here returns similar val- 
ues for the T] parameter at the level of few tens of percent 
at most (as determined by the nuclear radio fluxes), without 
violating the ferm/-LAT upper limits (see column 8 in Ta- 
ble 2 and the SEDs presented in Figures 8). Figure 4 presents 
the distribution of the mixing parameter 77 emerging from the 
SED matching. This distribution suggests that, even though 
the particular frequency ranges may be dominated by one of 
the two main emission components, the total observed lumi- 
nosities of the nuclear jets in BLRGs constitute on average 
not less than 1% of the accretion-related luminosities (me- 
dian 77 w 0.10), and that in the case of a few particular objects 
the jet-to-disk luminosity ratio may even approach unity. An 
interesting implication of the above statement, strengthened 
thanks to the inclusion of the ferm/-LAT data in our model- 
ing, is that one should expect a non-negligible non-thermal 
emission component in BLRGs in the MeV energy range, 
constituting as much as'-^ 1% — 10%of the total power 
emitted in the hard X-ray domain. If correct, this may be 
of importance for understanding the recently debated origin 
of the extragalactic MeV background (see the discussions in 
llnoue et al.lf2008t lAjello et al...2009) . On the other hand, the 
contribution of the jet emission to the total radiative output 
of the Seyfert 1 galaxies are in general very small (median 
77 w 5 X 10^'*; see Figure 4 and the SEDs shown in Figure 9). 

To investigate further the collected dataset, we plot the op- 
tical [i^ivle versus X-ray [7^F;y]2-iokcV fluxes for all the 
sources included in the sample in Figure 5 (the top panel), 
with BLRGs denoted by blue filled circles and Seyferts by 
red crosses. The bottom panel of Figure 4 presents instead the 
radio [t^Fi^Jsghz versus X-ray fluxes. Here, BLRGs are plot- 
ted as blue filled cycles when total radio fluxes are used, and 



Broad Line Radio Galaxies Observed with Fermi-LPd 



as cyan open circles when nuclear radio fluxes (obtained pre- 
dominantly from VLBI observations where available, supple- 
mented by VLA measurements in a few cases) are considered. 
In both panels, large symbols indicate the two objects detected 
by ferm/-LAT (3C 111 and 3C 120), while the medium-size 
symbols represent Pictor A galaxy which, even though not 
formally detected, is characterized by the third-highest TS in 
the analyzed two-year fermZ-LAT dataset. These flux-flux 
plots show that the optical and X-ray fluxes of BLRGs are 
roughly linearly correlated, as expected, though with a sub- 
stantial (order-of magnitude) scatter. In addition, however, it 
is revealed that 3C 1 1 1, 3C 120, and also Pictor A, are at the 
same time the brightest in radio, but only when the nuclear 
(and not the total) radio fluxes are taken into account. More- 
over, the nuclear radio fluxes seem well correlated with the 
X-ray fluxes for the whole BLRG population. Hence one may 
conclude that Ferm/-LAT detects preferentially those BLRGs 
which are characterized by the brightest radio and X-ray nu- 
clei. This supports the idea stated previously that the GeV 
emission of BLRGs is dominated by the innermost parts of 
their jets, and is therefore 'blazar-like,' being dependent on 
the jet luminosity and viewing angle. As argued above, the 
quasar 3C 273, for which the GeV luminosity is about 100 
times higher than the nuclear radio luminosity, should provide 
a reasonably accurate template for such an emission. Indeed, 
as illustrated in Figure 6, the 7-ray-to-radio energy flux ra- 
tios for the two BLRGs detected by Ferm/-LAT are of the or- 
der of ~ 100, while the corresponding upper limits for all the 
other objects from the sample (including Seyferts) are above 
or much above this value. This indicates that the detections of 
BLRGs in 7 rays are at present limited by the sensitivity of the 
ferm/-LAT instrument. In fact, from visual inspection of Fig- 
ure 8 and 9, the predicted 7-ray flux from the template is very 
close to the ferm/-LAT upper limits for Pictor A, 3C 390.3, 
3C 445, B3 0309H-41 IB and PKS 2153-69. 

Also it is important to note that the pc-scale components of 
the radio jets in both 3C 111 and 3C 120 are characterized 
by apparent superluminal motions with maximum velocities 
[3app :^ 5.9 and 5.3, res pectively, from MOJAVE monitoring 
data (iLister et al.ll2009h . The corresponding upper limits to 
the jet viewing angles are thus 9j < 20°. Interestingly, mildly 
relativistic velocities are inferred for 3C 3 90.3 and Pictor A 
(/3aDD up to c± 2.2 and 1.6, respectively; iKellermann et al.l 
120041; iTingav et al.ll2000l) . both of which are good candidates 
for future Fermi-LAT detections (see Table 2 and Figure 8). 
These four BLRGs, not coincidentally, have the brightest ra- 
dio/VLBI cores in the analyzed sample, consistent also with 
a relatively large degree of relativistic beaming (see the com- 
ment in AppendixA). Note also that 3C 303 and 3C 382 are 
othe r BLRGs with possible mildly relativistic pc-scale VLBI 
jets ( Giovannini et al.ll200l1) . 

But what else — besides the nuclear radio jets — could 
be the source of high-energy emission in the considered ob- 
jects? Clearly, large-scale structures (i.e., extended lobes, 
large-scale jets and terminal hotspots), which are particularly 
prominent in BLRGs, are expected to contribute at least at 
some level to the 7-ray emission in the GeV band, since the 
non-thermal synchrotron continua of these structures extend 
to the highest radio frequencies and even into the X-ray en- 
ergies in some cases (see in this context ferm/-LAT detec- 
tion of the e xtended lobes in t he low-power radio galaxy 
Centaurus A; lAbdo et ani2010el) . Such structures are on the 
other hand almost absent in Seyfert galaxies, due to much 
less prominent jet activity in those objects. Other possible 



7-ray emission sites in both BLRGs and Seyferts are their ac- 
cretion disks and disk coronae. Several related (and rather 
preliminary) studies presented in the literature, even though 
not directly applicable to the types of astrophysical objects 
discussed here, indicate that other additional processes than 
those typically considered for generating high-energy photons 
within the accreting matter (e.g., proton-proton interactions), 
may be relevant in this context. In these scenarios, the result- 
ing disk-related 7-ray emission may strongly depend on the 
main parameters of a black hole/accretion disk system (such 
as the black hole mass, sp in, accretion rate, and disk inclina- 
tion; see, e.g.,lMa hadevan et al.ll997trOka & Manmotol2003l: 
iNiedzwiecki et al...2009.) . Hence it is interesting to look into 
this issue in more detail for all the sources included in our 
sample. 

Figure 7 presents therefore the dependence of the 7-ray lu- 
minosities (or upper limits for such), expressed in the Ed- 
dington units (i^/LEdd) or as a ratio of the GeV and X-ray 
monochromatic luminosities (Lj/Lx), on (i) the black hole 
mass A^BH (top left panel), (ii) the accretion ratio £acc/-^Edd 
determined from the model fitting (top right panel), (iii) the 
proxy of the radio loudness parameter Lr/Lb [bottom left 
panel), and (iv) the observed proxy of the accretion rate 
Ls/L^dd (bottom right panel). In the figure, BLRGs and 
Seyferts are denoted by blue and red symbols, respectively. 
Large filled circles represent the two BLRGs detected by 
Fermi-LAI (3C 111 and 3C 120), while the medium-size open 
circles represent Pictor A. In the bottom left panel, when nu- 
clear radio fluxes are used instead of the total radio fluxes, 
BLRGs are denoted by cyan symbols, and Seyferts by ma- 
genta symbols. A ratio of the GeV and X-ray monochromatic 
luminosities (L-f/Lx) expected from a hybrid model is given 
as dashed lines in the bottom right panel, for different values 
of a mixing parameter 77 = 10^^, 10^^, and 1. 

Note that the monochromatic luminosity ratio Ljx/Lb is 
simply proportional to the 'standard' radio-loudness param- 
eter TZ, and in particular that L^/Lb — {i^r/i^b) x 7?, '^ 
10^^ 7?.. Also, for the high accretion-rate objects analyzed 
in this paper one expects the accretion luminosity iacc — 
Ltot — 10 X Lb. As shown, the two galaxies which are the 
brightest in 7 rays, and also Pictor A, are not characterized 
by any outstanding values of the radio loudness, accretion 
rate, or black hole mass, but are in fact quite representative 
for the other BLRGs included in the sample. This supports 
once again our conclusion that the detected GeV emission of 
BLRGs is related predominantly to nuclear jets, and not the 
other possible nuclear emission components. What should be 
noted for constraining theoretical models regarding the high- 
energy emission of accretion disks in AGN is the fact that, 
at least for a number of sources, the Fermi-LAiT upper limits 
evaluated here already probe the GeV emission of BLRGs and 
Seyferts down to the levels of 1% of the X-ray (disk-related) 
luminosity, or equivalently 0.01% Eddington luminosity (see 
Figure 7). 
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APPENDIX 

Here we present all the SEDs of the BLRGs (Figure 8) and Seyfert 1 galaxies (Figure 9) analyzed in this paper which are not 
detected by Ferm/-LAT , together with the hybrid model fits described in the main text. The figures illustrate the expected level 
of the Ge V emission of each source considered, in comparison with the current upper limits provided by two-years of Fermi-LAT 
data. Note in particular that the expected 7-ray fluxes of Pictor A, 3C 390.3 3C 445, B3 0309H-41 IB and PKS 2153-69 are close 
to the Ferm/-LAT limits derived in this paper (cf.. Table 2). In contrast to the BLRGs (or at least the brightest and most beamed 
examples of such considered in this work), we speculate that the jet-related non-thermal emission of all the analyzed Seyfert 1 
galaxies within the GeV photon energy range are in general beyond the level of detectability with Fenni-LAI . 
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TABLE 1 
MULTIWAVELENGTH PROPERTIES OF THE ANALYZED BLRGS AND SEYFERTS 
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13,22,25,54,61-63,76.77 

13,22,25,54,64,65,76 

13,22,25,54,64-66,76 

9,13,22,25,54, 

61,64,67,76 

13,22,25,54,68-71,76 

9,22,25,54,61,64,72,76 

9,22,25,54,61,68,71,76 

9,25,54,61,64,74,76 



^ References: [1 1 Gnindi et al. : 2006); |21 Hewitt & Biii-bidge 1 1991 ); |3| Laiirent-Muehleisen et af ( 1997); |4| Marchesini et al. : 2004); |51 Brinkmann et al. : 1994); [6] Molina et al. ( 2008); |71 Molina et al. 1 2009); [8] Henstock et al. ( 1995); 
[9]IBirdetal. (2007); [10] Wozniak et al. 1 1998); [11] Sarabriina et al. (1999); |12] Eracleous et al. (2000); [13] Tiieller et al. (2010); [14] Turner & Poiind,s (1989); [15] Linfield (1987); [16] Lewis et al. (2005); [17] Gliozzi et al. (2003); 
[18] Kataokaetal. (2007); [19] Hirabavashi et al. (2000); [20] Fomalont et al. (2000); [21] Zhou et al. (2010); [22] Fukazawa et al. (2011); [23] Ogle et al. (2005); [24] Ballantyne et al. (2004); [25] Sikora et al. (2007); [26] Tingay et al] 
(2003); [27] Tingay et al. ( 2000); [28] Eracleous & Halpern (1998); [29] Hough et al. ( 2002); [30] Labita et al. ( 2006); [31] Reeves & Turner ( 2000); [32| Crawford & Fabian ( 1995); [33] Hardcastle et al. ( 2007); [34] Ajello et al. ( 2008); [35] 
[feataoka et al. ( 2003); [36] Gioyannini et al. ( 1990); [37] Donato et al. ( 2005); [38] Gioyannini et al. ( 1994); [39| Sambruna et al. ( 201 1); [40] Gliozzi et al. ( 2009); [41] Sambruna et al. ( 2009); [42] Pearson et al. ( 1992); [43] Wang et al. ( 2009); 
[44]Tayloretal. (1996); [45] Tazaki et al. ( 2010); [46| Tingay et al. (2002); [47] Young et al. (2005); [48] Preuss & Fosbui-y ( 1983); [49] Grandi et al. ( 2007); [50] Sambruna et al. ( 2007); [511 Sambruna et al. ( 1998); [52] Kellermann et al] 
[1989); [53] Liu et al. (2006); [54] Ka,SDi et al. (2005); [55| EdeLson (1987); [56] Gallo (2006); [57] Inoue et al. (2007); [58] Larsson et al. (2008); [59| Bianchi et al. (2001); [60| O'Neill etal. (2007); [61| Shinozaki et al. (2006); [62] 
[SchmoUetal. (2009); [63] Gondoin et al. (2001); [64] Nandra et al. (2007); [65] Sazonoy et al. ( 2007); [66] Dasgunta & Rao (2006); [67] De Rosa et al. ( 2002); [66] Beckmann et al. (2006); 1691 Rey nolds et ak a004.); l701 .Brenneman et al] 
tSOOTI) ; [711 IMarkowitz & Reevej|2009l) ; [721 ISteenbrugge et alj 120051) ; [731 IDe Rosa et al J I2004[) ; [741 IBlustin et alj 120031) ; [751 INeff et alj|1995l) ; [761 IHcl 120021) ; [771 IVeron-Cettv et alj|1991h . 

Flux estimated from the counts rate of Chandra ACIS CCD chip using PIMMS and assuming the X-ray photon index Fx — 1.5. 
'^ Unreliable extinction estimate due to the location of the source at low Galactic latitude | fc | < 5 deg. 



5 GHz flux was obtained from our new VLBA obseryations obtained on 2009 Feb 7 in program BC 185 (PI; C. C. Cheung). 
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TABLE 2 
Results of the Fermi-hAT data analysis 



-P'>0.1GcV 
[10-'-'phcm-2s-i] 



[»^-Fi/]o.l-10GoV 



logL^ 
[ergs^i] 



log Lace" 

[ergs^i] 



BLRGs 



3C18 


1.5 


2.5^^ 


<9.2 


B3 0309+41 IB 


< 1 


2.5f 


<8.9 


3C 111 


31 


2.7±0.2 


35±12 


3C 120 


34 


3.0±0.3 


37±14 


Pictor A 


20 


2.5f 


<15 


3C215 


13 


2.5f 


< 14 


3C227 


< 1 


2.5f 


<5.0 


3C 303 


2.0 


2.5f 


<5.9 


RGB J 1722+246 


< 1 


2.5f 


<38 


3C382 


2.8 


2.5f 


<12 


3C 390.3 


3.2 


2.5f 


<7.4 


3C411 


8.4 


2.5f 


< 18 


4C 74.26 


3.1 


2.5f 


<11 


S5 2116+81 


1.6 


2.5f 


<5.9 


4C 50.55 


< 1 


2.5f 


< 180 


PKS 2153-69 


1.4d 


2.5f 


<10 


3C445 


< 1 


2.5f 


<2.1 


PKS 2251+11 


< 1 


2.5f 


<5.4 



<4.0 
<3.9 

15 

16 
<6.3 
<6.1 
<2.2 
<2.6 
< 16 
<5.2 
<3.2 
<7.6 
<4.7 
<2.6 
<77 
<4.4 
< 0.90 
<2.3 



C44.6 


45.3 


C44.3 


45.5 


43.9 


45.0 


43.6 


44.9 


C43.2 


44.5 


C45.6 


45.9 


C43.6 


44.9 


C44.1 


45.0 


C45.1 


45.0 


C43.6 


45.1 


C43.4 


45.0 


C45.8 


45.7 


C44.1 


45.6 


C43.6 


45.3 


C43.8 


44.4 


C42.9 


44.2 


C42.8 


44.8 


C44.9 


45.4 



0.070 

0.12 

0.35 

0.15 

0.25 

0.034 

0.017 

0.20 

0.058 

0.025 

0.075 

0.17 

0.015 

0.023 

0.024 

0.10 

0.022 

0.010 



Seyfert.s 



2.6x10-"' 
< 4.0x10-4 
2.3x10-* 
3.3x10-4 
4.6x10-* 
1.3x10-* 
1.0x10-3 
3.2x10-* 
1.2x10-3 



Mrk 335 


8.9 


2.5f 


Fairall 9 


<1 


2.5f 


Ark 120 


< 1 


2.5f 


Mrk 110 


3.9 


2.5f 


NGC 3783 


< 1 


2.5f 


NGC 4593 


2.5 


2.5f 


IC4329A 


15 


2.5f 


Mrk 509 


< 1 


2.5f 


NGC 7469 


< 1 


2.5f 



< 11 

< 1.9 

<2.8 
<6.9 
<3.2 

< 8.9 

< 19 

<2.7 
<7.1 



<4.8 
<0.81 

< 1.2 
<3.0 

< 1.4 
<3.9 
<8.1 

< 1.2 
<3.1 



<42.9 
<42.6 
<42.5 
<42.9 
<41.5 
<41.8 
<42.7 
<42.5 
<42.3 



44.6 
45.0 
44.9 
44.8 
44.1 
43.9 
44.7 
45.1 
44.8 



^ Accretion luminosity derived from the SED fitting, assuming hybrid model discussed in the paper. 

Mixing parameter described in the text, derived from the SED fitting. 
^ Photon spectral index was fixed at 2.5. 

A nearby LAT source was found 0.27 deg apart from PKS 2153-69 with TS of 44. We treat this as a background source since its separation from the target is more than 2 times larger 



than its 95 % 7-ray localization error radius. 
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Fig. 8. — Broad-band SEDs of the BLRGs which are not detected at high significance in the GeV photon energy range. Fermi-hAT upper Hmits are indicated 
by red arrows. Black squares represent the historical data from NED. Magenta squares denote the 5 GHz radio fluxes of the unresolved nuclei. The green curves 
con'espond to the template of the accretion-related Seyfert-type emission (from Koratkar & Blaes 1999), matched to the infrared-to-X-ray continuum of each 
source. The blue curves coiTespond to the broad-band spectrum of the quasar 3C 273 (from Soldi et al. 2008), used here as a template of the jet-related emission 
and scaled to match the radio fluxes for each source. The mixing parameter 77 for the phenomenological hybrid model discussed in § 4 is given in each panel. 
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Fig. 8. continued. 
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Fig. 8. continued. 



Broad Line Radio Galaxies Observed with Fermi-LAJ' 



17 



iJ' in-14- 



iJ' in-14- 



-^ . 



Mrk 335 
(7)=2.6x10"*) 




10-° 10"* 10" 10' 10^ 10'2 

Energy [eV] 



Ark 120 
(7)=2.3x10"*) 




iJ" 10-1* 



10" 10' 10° 10' 

Energy [eV] 



NGC 3783 
(77=4.6x10"*) 




uj" lo-l* 



Energy [eV] 



*. . 



Fairall 9 
(7)<4.0xl0"*) 



J 




10-8 ,0-* 10° 10* 10^ 10'2 

Energy [eV] 



Mrk 110 
(77=3.3x10"*) 




10" 10' 10° 10' 

Energy [eV] 



NGC 4593 
(77=1.3x10"') 




Energy [eV] 



Fig. 9. — Broad-band SEDs of the high-accretion Seyfert 1 galaxies. Fermi-hAT upper limits are indicated by red aiTows. Black squares represent the histoiical 
data from NED. Magenta squares denote the 5 GHz radio fluxes of the unresolved nuclei, except for Fairall 9 where a radio detection has not been reported in the 
literature. The green curves correspond to the template of the accretion-related Seyfert-type emission (from Koratkar & Blaes 1999), matched to the infrared-to- 
X-ray continuum of each source. The blue curves coiTespond to the broad-band spectrum of the quasar 3C 273 (from Soldi et al. 2008), used here as a template 
of the jet-related emission and scaled to match the radio fluxes for each source. The mixing parameter t] for the phenomenological hybrid model discussed in § 4 
is given in each panel. 
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Fig. 9. continued. 



